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SATELLITE-BASED POSI TIONING SYSTEM IMPROVE MENT 

This application claims the priority filing date of U.S. Provisional Application 
Serial No. 60/416,367 filed on October 4, 2002. 

Field of the Invention 

This Invention relates to the design of receivers employed in satellite-based 
positioning systems (SPS) such as the US Navstar Global Positioning System 
(GPS), the Russian Global Navigation Satellite System (GLONASS) and the 
European Galileo system. More specifically, the invention relates to methods, 
devices and systems for determining a receiver location using weak signal satellite 
transmissions. 

Background of the Invention 

Satellite based positioning systems operate by utilizing constellations of 
satellites which transmit to earth continuous direct sequence spread spectrum 
signals. Receivers within receiving range of these satellites intercept these signals 
which carry data (navigation messages) modulated onto a spread spectrum carrier. 
This data provides the precise time of transmission at certain instants in the signal 
along with orbital parameters (e.g., precise ephemeris data and less precise 
almanac data in the case of GPS) for the satellites themselves. By estimating the 
time of flight of the signal from each of four satellites to the receiver and computing 
the position of the satellites at the times of transmission corresponding to the 
estimated times of flight It is possible to determine the precise location of the 
receiver's antenna. 

In a conventional SPS receiver, the process by which this is done involves 
estimating pseudoranges of at least 4 satellites and then computing from these the 
precise location and clock error of the receiver Each pseudorange is computed as 
the time of flight from one satellite to the receiver multiplied by the speed of light and 
is thus an estimate of the distance or 'range' between the satellite and the receiver. 
The time of flight Is estimated as the difference between the time of transmission 
determined from the navigation message and the time of receipt as determined using 
a clock in the receiver. Since the receiver's clock will inevitably have a different 
present time when compared to the clock of the satellites, the four range 
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computations will have a common error. The common error is the error in the 
receiver's clock multiplied by the speed of light. 

By using at least 4 satellites it is possible to solve a set of equations to 
determine both the receiver clock error and the location of the antenna. If only 3 
5 measurements are available it is still possible to determine the location and clock 

error provided at least one of the receiver's coordinates is already known. Often, this 
situation can be approximated by estimating the altitude of the antenna. 

The signals from the satellites consist of a carrier signal which is biphase 
modulated by a pseudo-random binary spreading code at a relatively high "chipping" 

10 rate (e.g.. 1 .023 MHz) and then biphase modulated by the binary navigation 

message at a low data rate (e.g., 50 Hz). The carrier to noise ratio is typically very 
low (e.g., 31dBHz to 51dBHz) at the earth's surface for a receiver with unobstructed 
line of sight to the satellite from its antenna. However, it is sufficient to permit the 
signals to be detected, acquired and tracked using conventional phase-locked loop 

15 and delay-locked loop techniques and for the data to be extracted. 

The process of tracking the code of a signal in a conventional SPS receiver 
involves the use of a hardware code generator and signal mixer. When the locally 
generated code is exactly aligned with that of the incoming signal, the output from 
the mixer contains no code modulation at all. Hence the bandwidth of the signal is 

20 much less and it can be filtered to greatly increase the signal to noise ratio. This is 
usually done using a decimation filter such that the correlator output sampling rate is 
much lower than the input sampling rate (e.g., 1kHz at the output compared to 
1.3MHz at the input). 

Also, in the case of GPS, the precise time of transmission of this signal 

25 corresponding to any given instant at the receiver can be determined by latching the 
state of the code generator to get the code phase and by counting the code epochs 
within each bit of the data and by counting the bits within each word of the navigation 
message and by counting the words within each subframe of the message and by 
extracting and decoding the times of transmission corresponding to the subframe 

30 boundaries. A similar scheme can be used for any SPS. 

However, traditional SPS receivers can suffer from troublesome lapses in 
position identification in the presence of weakened transmission signals. When the 
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direct line of sight between the antenna and the satellites is obstructed, signals may 
be severely attenuated when they reach the antenna. Conventional techniques can 
not be used to detect, acquire and track these signals. Moreover, under these 
circumstances even if the signal could be detected, the carrier-to-noise ratio of a 
5 GPS signal, for example, may be as low as or lower than 24dBHz and as such it is 
not possible to extract the data from the signals. 

Prior art devices have attempted to minimize or overcome these shortcomings 
through the use of aiding information. In such schemes, additional information is 
externally supplied to the SPS receivers through various secondary transmission 

10 sources to balance the shortfall of information resulting from the attenuated signals. 
Examples of such devices are taught in the patents to Taylor et al. (US Patent No. 
4.445,1 18) (aided by satellite almanac data); Lau (U.S. Patent No. 5,418,538) (aided 
by differential satellite positioning information and ephemerides); Krasner (US Patent 
No. 5,663,734) (aided by transmission of Doppler frequency shifts); Krasner (US 

15 Patent No. 5.781 ,156) (aided by transmission of Doppler frequency shifts); Krasner 
(US Patent No. 5,874.914) (aided by Doppler, initialization and pseudorange data) 
Krasner (US Patent No. 5,841,396) (aided by satellite almanac data); Loomis, et al. 
(US Patent No. 5,917,444) (aided by selected satellite ephemerides, almanac, 
ionosphere, time, pseudorange corrections, satellite index and/or code phase 

20 attributes); Krasner (US Patent No. 5,945,944) (aided by timing data); Krasner (U.S. 
Patent No. 6,016,1 19) (aided by retransmission of data from satellite signal) 

However, aiding information requires additional transmission capabilities. For 
example, aiding information may be sent to the SPS receiver using additional 
satellite transmitters or wireless telephone systems. As such, it is a significant 

25 advantage to reduce the quantum of aiding information supplied to limit the use of 
such additional resources. For example, when the voice path of a wireless 
communication network is being used to communicate the aiding information, the 
voice communication will be interrupted by the aiding message. The aiding 
messages must therefore be as short as possible in order to limit the voice 

30 interruptions to tolerable durations and frequencies. Also, no matter how the aiding 
data is communicated, its communication will delay the operation of the receiver. In 
many applications the location data is needed promptly and therefore any delay must 
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be minimized. 

Brief Description of the Invention 

The present invention is an improvement on the invention disclosed in 
5 International Patent Application PCT/AU01/00519. 
Calibration of Locai Osciiiator 

\ in many assisted GPS applications rapid acquisition of satellite signals is a 
key requirement. Acquisition is delayed because of drift of the reference oscillator in 
the receiving unit. The relative velocity along the line of sight from the receiver to the 

10 satellite induces a Doppler shift in the frequency of received signal. The Doppler 
shift contains useful information on the velocity of the receiver antenna, but the 
presence of the Doppler shift necessitates a frequency search that increases the 
time for acquisition. Reference oscillator drift is a major contributor to lengthening 
acquisition time as it causes the "Doppler" frequency search to be increased to allow 

15 for reference oscillator drift. By utilizing the precise signal framing of a digital 

communications link, the invention calibrates a local oscillator and thus reduce the 
effect of drift. This is accomplished by counting local oscillator cycles and fractions 
thereof over a period precisely determined by a number of signal framing intervals. 
Once the calibration offset is determined it is used as a correction by the GPS 

20 receiver firmware when performing acquisition searches or it can be used to correct 
the oscillator frequency so as to minimize the offset. 

Cancellation of Cross Correlation 

Another aspect of the Invention is the reduction of cross correlations between 
25 weak and strong signals experienced at correlator outputs. These cross- 
correlations are inherent limitations of the GPS C/A Code structure. The cross 
correlations between codes at certain code and Doppler offsets are only 20dB 
between the peak of the autocorrelation main lobe. At the correlator output these 
cross-correlations are indistinguishable from correlations with the locally generated 
30 replica of the weak signal being sought. Reducing the level of cross-correlations 
caused by a strong signal will reduce its jamming effect on weaker signals. In this 
way the usable dynamic range is increased to permit weak signals to be acquired, 
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tracked and used in the presence of strong ones. At least 3 satellites are needed to 
make a 2D fix and 4 satellites are require to make a 3D fix. Furthermore, more than 
the minimum number may be required to obtain a low enough Dilution Of Precision 
to permit an accurate fix to be made. Hence the ability to use more of the signals 
present is an advantage. In urban canyons this advantage will be distinct in that 
there are often only one or two strong signals present and these jam all of the 
weaker ones. The biggest problem with the concept of canceling the strong signals 
is that the signals are represented with very low precision at the input to a correlator 
and, hence, any scaling of the signal can only be extremely crudely performed. This 
threatens the viability of the concept. The present invention provides allows the 
scaling to be performed at a point where the signal is represented by 10 bit samples 
and scaling is much more feasible. 
Doppter Location Scheme 

Another aspect of the invention deals with the problem that in a weak signal 
environment, it is not possible to resolve code phase ambiguity without knowing the 
GPS receiver's initial position to within 150km. The 150 km is the distance a satellite 
signal would travel in the time occupied by one half of a code epoch. The invention 
provides a method for computing the initial position autonomously without requiring 
prior knowledge of the location by computing it from measured satellite Doppler 
differences. Using differences to perform the calculation ensures that any 
dependence on the current local oscillator offset is removed from the calculation. 
Ensemble Averaging Scheme 

A still further aspect of the invention concerns the limitations on sensitivity 
caused by limitation in performing Fourier transforms. The transform involves 
accumulating values in frequency bins. Previous algorithms are limited in sensitivity 
because the reduction of binwidth is offset by squaring losses. This causes growth 
of the FFT (fast Fourier transform) size to impractical lengths. To overcome this 
problem a spectral averaging algorithm is employed for integration periods beyond 
160ms (FFT length of 128), To avoid excessive memory usage for storing FFT 
arrays, the invention operates in the fashion of a moving average by filtering the 
squared magnitudes of the FFT bins in the same fashion as for the autocbnvolutions 
when acquiring. 
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Calculation of Almanac Aiding from Table of Orbit Coefficients 

A still further aspect of the invention reduces the traffic from the receiver to 
the aiding source. In particular, an alternative is provided to having the receiver 
acquire its own almanac data and requesting augmentation coefficients for a specific 
5 Issue Of Almanac as aiding information. Instead the invention hardcodes orbital 

coefTicients into a lookup table and the aiding information is provided as broadcasts 

at regular intervals. 

Absolute Time Determination 

Another aspect of the invention concerns the problem of extracting data from 

10 the GPS signal when there is a poor signal to noise ratio. Often the output of a weak 
signal GPS receiver is limited to code phases for each satellite rather than a full 
time-of-transmit (pseudorange). However a full time-oMransmit is obtainable from a 
partial code phase provided that a reference time stamp for the observation is 
available, as well as an estimate of the users location. The estimate of the users 

15 location can be used to calculate a nominal time of flight for each satellite provided 
that its uncertainty is less than 150 km (0.5 ms), while the reference time stamp is 
used to convert the time-of-receipt into an estimated satellite time-of-transmit. The 
code phase is then used to refine the time-of-transmit estimate into an ambiguity 
resolved time-of-transmit. This ambiguity resolved time-of-transmit will be consistent 

20 with the true user position and hence can be used to calculate the true user position. 
Modified Search Engine 

A still further aspect of the invention is a modified search engine. Modem 
GPS receivers often incorporate massively parallel correlation hardware to speed the 
acquisition process by concurrently searching across a broad range of codephases 

25 and/or Doppler frequencies. Such hardware are sometimes referred to as search 
engines. The signal processing algorithms utilitzed in the invention can be 
incorporated into a modified search engine to achieve even faster searching by 
reducing the length of integration period required. 

An object of the present invention is to provide a satellite-based positioning 

30 system receiver for weak signal operation have each of the advantages just 
described. 

Consistent with these objectives, a device made in accordance with this 
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invention utilizes a novel signal processing scheme for detecting, acquiring and 
tracking attenuated satellite signals, such as those that might be received at an 
indoor location, and computes location solutions. The scheme makes novel use of 
attenuated satellite signals and minimal externally-supplied aiding information. 

5 

Brief Description of the Drawings 

This invention is illustrated by means of the accompanying drawings. 
However, these figures represent examples of the invention and do not serve to limit 
its applicability. 

10 FIG. 1 is a sequence diagram describing the interactions between an aiding 

source, a call taker and a handset with an integrated SPS receiver according to one 
embodiment of this invention; 

FIG. 2 is a flowchart describing the overall algorithm according to one 
embodiment of this invention for acquiring satellite signals, measuring code phases, 
15 carrier smoothing these measurements, computing pseudorange differences and 
computing handset location; 

FIG. 3 is a block diagram of a typical SPS receiver according to this 
invention.; 

FIG. 4 is a block diagram describing the signal processing algorithm used to 
20 measure amplitude in each of an early and a late arm of each channel of the 
correlator according to one embodiment of this invention; 

FIG. 5 is a block diagram describing the carrier smoothing algorithm used to 
reduce the error in the code phase measurements according to one embodiment of 
this invention; 

25 FIG. 6 is a block diagram describing the algorithm used to compute the 

location and velocity from the code phase and carrier frequency differences 
according to one embodiment of this invention. 

FIG. 7 is a schematic diagram of a canceler Integrated into a correlator. 
FIG. 8 illustrates a 2 bin search engine architecture. 
30 FIG. 9 is a schematic representation of a 2 codephase modified search 

engine. 
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Detailed Description of Embodiments of the Invention 

This invention relates to refinements and extensions to a commonly owned 
invention disclosed in U.S. Patent No. 5,459,473. Accordingly, the foregoing U.S. 
patent is hereby incorporated by reference. 
5 Aiding Source/Receiver Interaction 

As previously described, the aiding data used in accordance with the present 
invention may be limited to information that includes an approximate location for an 
SPS receiver and the positions and velocities of a specific set of satellites. This 
information is determined and provided through a request/response sequence. A 

10 model of one embodiment of such an exchange in accordance with the present 
invention is depicted in FIG. 1 . 

A typical exchange might involve an SPS Receiver 1 , an Aiding Source 2 and 
a Call Taker 3. For Instance, the SPS Receiver 1 might be a GPS receiver 
embedded in or co-located with a wireless telephone or other handset. The Aiding 

15 Source 2 may be located at a call center or cell site or elsewhere in the wireless 

network such that the aiding data is transmitted via a wireless communication link to 
the handset. The Call Taker 3 may also be located at the call center or other 
location accessible from the wireless network. The ultimate user of the location data 
may be either the Call Taker 3 or the user accompanying the SPS Receiver 1 . Other 

20 forms of transmission between the SPS Receiver 1 , Aiding Source 2, and the Call 
Taker 3 may be utilized without departing from the objectives of the present 
invention. 

To begin the exchange, the SPS Receiver 1 sends a First Aiding Request 4 to 
the Aiding Source 2. This would typically occur upon activation of the SPS Receiver 

25 1 but may occur at other times as well In response, the Aiding Source 2, sends a 
First Aiding Response 5 which contains the approximate location of the SPS 
Receiver 1 . Preferably, the approximate location of the SPS Receiver 1 is accurate 
to better than one half of a code epoch of a satellite signal multiplied by the speed of 
light or about 150km in the case of GPS. The approximate location may also be sent 

30 to the Call Taker 3 in a First Aiding Report 6. 

DoPDier Location Scheme Alternative 
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In a weak signal environment, it is not possible to resolve code phase 
ambiguity without knowing the GPS receiver's initial position to within 150km. A 
method for computing the initial position autonomously without requiring prior 
knowledge of the location is to compute it from measured satellite Ooppler 
5 differences. Using differences to perform the calculation ensures that any 

dependence on the current local oscillator offset is removed from the calculation. A 
Doppler positioning algorithm that implements this scheme is the following, which 
would be executed by microprocessor 20 of Figure 3: 

1 . For each input observation made by the GPS receiver apply the space vehicle 
10 (SV) clock corrections to both the satellite time of transmits (Toti*) and the satellite 

Doppler frequencies. 

2. For each satellite 'i' calculate satellite position XI (Xi.Yi,Zi) at the corrected 
satellite time of transmit Toti'. 

3. Calculate the average satellite position by taking the vector mean of the 
15 calculated satellite positions. Convert the vector mean to a unit vector and then 

scale by the radius of the earth. This provides an initial estimate for the user position 
Xu (Xu,Yu,Zu). The user's local velocity is assumed to be zero throughout the entire 
process, so (XuDot,YuDot,ZuDot) is (0,0,0). 

4. Iterate through the following steps 5 to 10 until the residual error is sufficiently 
20 small (say less than 100 m). 

5. Using the current estimate of the user position and user velocity (zero), 
calculate the range 'Ri' and range-rate 'RRf vectors to the satellite 'i'. Note that 
these will change at each iteration as the estimate of the user position Xu converges 
to the true position. 

25 6. Using the range vector 'Ri' calculate the line-of-sight vector 'LOSr (the unit 
range-vector). Calculate the projection of the range-rate vector onto the line of sight 
vector (RRi . LOSi) thereby determining an estimate of the predicted satellite range- 
rate. Finally calculate a range-rate prediction vector 'yh as the vector of differences 
between the first satellite range rate and all remaining range-rates. 

30 7. Using the measured satellite observations, calculate a measured range-rate 
vector y as the differences between the first measured satellite range-rate (Doppler 
frequency scaled to same units as predictions above). 
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8. Using the satellite line-of-sight vectors calculate a measurement matrix 'M' of 
line of sight vector differences between the first satellite and all remaining satellites. 
Each row in the matrix corresponds to the difference between the first satellite LOS 
vector and the one of the remaining satellite LOS vectors. 
5 9. Calculate the pseudo-inverse of the M matrix 'Mpi'. 

10. Correct the current estimate of the user position by calculating dx = Mpi*(y- 
yh) and Xu = Xy + dx. When the magnitude of the correction *dx' is less than a 
sufficiently small threshold (say less than 100 m) then terminate the iteration 
process. 

10 

Absolute Time Determination 

Related to the Doppler position process is the following aspect of the 
invention. One of the limitations of assisted GPS in a weak signal environment is the 
inability to extract data from the GPS signal due to the poor signal to noise ratio. 
15 This means that the output of a weak signal GPS receiver is often limited to code 

phases for each satellite rather than a full time-of-transmit (pseudorange). However 
a full time-of-transmit can be obtained from a partial code phase provided that a 
reference time stamp for the observation is available, as well as an estimate of the 
users location. 

20 The estimate of the users location can be used to calculate a nominal time of 

flight for each satellite provided that its uncertainty is less than 150 km (0.5 ms), 
while the reference time stamp is used to convert the time-of-receipt into an 
estimated satellite time-of-transmit. The code phase is then used to refine the time- 
of-transmit estimate into an ambiguity resolved time-of-transmit. This ambiguity 

25 resolved time-of-transmit will be consistent with the true user position and hence can 
be used to calculate the true user position. 

For the above process to provide acceptable positional accuracy it Is 
necessary that the reference time-stamp have an error of less than 10 ms. The 
algorithm that is provided below provides a means by which any error in the 

30 estimated time-of*receipt can be calculated providing that at least one additional 
satellite is available for use. This algorithm is similar in concept to the Doppler 
position process except that full or ambiguity resolved time*of-transmits are 
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employed. 

The algorithm requires an additional redundant satellite because attempting to 
estimate the time-of-receipt error is equivalent to introducing an additional degree of 
freedom or an additional unknown and hence an additional linearly independent 
observation is required for this to be estimated. If only 4 obsen/ations are available 
then there will only be a single position solution to the problem regardless of the 
time*of-receipt and the time-of-receipt error will always be calculated as zero. 

The algorithm is the following: 

1. For each raw input observation from satellite i, correct the satellite time-of- 
transmit Toti for space vehicle (SV) clock corrections. 

2. For each satellite I calculate the satellite position vector Xi (Xi.Yi.Zi) and 
satellite velocity vector XiDot (XiDot,YiDot.ZiDot) at the corrected Toti. The 
relativistic corrections for each satellite can also be calculated at this stage. 

3. Apply relativistic corrections to the corrected Tot's and recalculate SV position 
and velocity vectors. 

4. V- Calculate a starting point for the user position Xu (Xu.Yu.Zu) by taking the 
mean satellite position, converting it to a unit vector and scaling it up by the radius of 
the earth. 

5. Iterate through the following steps 6 to 19 until the residual position error is 
sufficiently small (say less than 100 m) and the residual time-of-receipt estimate 
error is also sufficiently small (say less than 10 ms). 

6. For each raw input observation from satellite i, correct the satellite time-of- 
transmit Toti for space vehicle (SV) clock corrections. Also apply relativistic 
corrections, (i.e. Steps (1) and (3)). 

7. Using the corrected Toti's for each satellite, recalculate the satellite position Xi 
and velocity XiDot vectors. 

8. Using the most recent user position Xu recalculate the satellite range Ri and 
range-rate RRi vectors. Corrections for earth rotation are also included within the 
range and range-rate vectors. 

9. An optional step at this point Is to is to apply additional corrections to the 
measured observations that take into account propagation delay effects caused by 
the ionosphere and troposphere. 
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10. Using the corrected satellite observations produce a vector yd of time-of- 
transmit (or pseudorange) differences scaled to units of meters. The differences are 
taken between the first Tot and each subsequent Tot. 

1 1 . Using the calculated satellite range vectors calculate a vector of range 

5 difference predictions (in units of meters) ydh for each satellite. The differences are 
taken between the first range prediction and each subsequent range prediction. 

12. Calculate a measurement (attitude) matrix M, where each row in the matrix is 
the difference between the first LOS vector and each subsequent LOS vector. 
Calculate the pseudo-inverse Mpi of the measurement matrix M. 

10 13. Calculate the correction to the user position Xu as dx = Mpi*(yd-ydh). 

14. If the number of algorithm iterations exceeds 3 iterations then include the 
following additional steps 15 to 19. 

15. Recalculate the satellite range Ri and range-rate RRi vectors using the new 
estimate for the user position Xu (calculated in step 13). 

15 16. Using the corrected satellite observations produce a vector y of satellite 

pseudoranges scaled to units of meters. This is similar to step (10) except that the 
vector y is a vector of full pseudoranges rather than pseudorange differences. Note 
that this calculation requires an estimate of the time-of-receipt which could contain 
an error that can be solved for. 

20 17. Using the calculated satellite range vectors calculate a vector of range 
predictions (in units of meters) yh for each satellite. This is similar to step (11) 
except that the vector yh is a vector of full pseudorange predictions. 

18. Calculate a 2"^ measurement matrix Ml consisting of 2 rows and Nobs 
columns, where Nobs is the number of satellite observations. The 1^* column 

25 contains the value c, where c is the speed of light in meters per second. The second 
column contains estimate of the satellite range-rate in units of meters per second, 
where the satellite range-rate is calculated as the projection of the satellite range- 
rate vector onto the line-of-sight vector. 

19. Calculate a dt vector as dt = pseudoinverse(M1)*(y-yh). The second element 
30 of dt is then an estimate of the time-of-receipt error given the value of the time-of- 
receipt used to calculate the pseudorange and pseudorange estimates. This time 
can then be used to correct the time-of-receipt of the raw observation and the 
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Further Aiding Source/Receiver Interaction 

With the approximate location and previously stored almanac data, the SPS 
5 Receiver 1 performs its correlation search to acquire satellite signals. The almanac 
data and the approximate location help to constrain the initial search once at least 
one satellite has been acquired. Upon acquisition, the SPS Receiver 1 sends a 
Second Aiding Request 7 to the Aiding Source 2. The Second Aiding Request 7 
includes infomiation for identifying the specific set of satellites used by the SPS 

10 Receiver 1 in determining pseudorange differences. In response, the Aiding Source 
2 determines the precise positions and velocities of the identified set of satellites 
from ephemeris data for the satellites. The determined positions and velocities are 
then sent to the SPS Receiver 1 in a Second Aiding Response 8. Since this elapsed 
time is known and assuming that the latency between transmission and reception of 

15 the request for aiding can be determined it is possible for the aiding source to 
determine the time of reception of the satellite signals to within a few tens of 
milliseconds. 

Therefore, under this scheme since the Aiding Source 2 provides precise 
satellite positions and velocities, the Aiding Source 2 rather than the SPS Receiver 1 

20 needs to be able to determine specific time synchronization data from the satellite 
signals and needs to maintain or acquire ephemeris data. Moreover, in order to 
ensure that the satellite positions are accurate when received by the SPS Receiver 1 
from Aiding Source 2, the latency period for the communication between the two 
must be within a few tens of milliseconds. This will ensure a limitation on the error in 

25 the computed satellite locations to a few meters. To this end, in the preferred 

embodiment, the Second Aiding Request 7 occurs at a known elapsed time from the 
instant when the code phases of the satellite signals latch. Since this elapsed time Is 
known and assuming that the latency between transmission and reception of the 
request for aiding can be determined it is possible for the Aiding Source 2 to 

30 determine the time of reception of the satellite signals to within a few tens of 
milliseconds. 

After receiving the satellite positions and velocities and using pseudorange 
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and range rate differences, the SPS Receiver 1 then computes a Position and 
Velocity (PV) solution to determine its precise location, speed, heading, etc. After 
such determination, the SPS Receiver 1 then sends a Receiver Report 9 to the 
Aiding Source 2 which includes the raw location, speed, heading, height, satellite 
5 identifications and the solution mode used by it (i.e. 3D or 2D with altitude aiding). 

Upon receipt of this later information, the Aiding Source 2 may take further 
action. For example, the Aiding Source 2 may use the known satellite set and times 
of transmission, select differential pseudorange corrections (obtained by any 
available means) and then compute a corresponding location correction consistent 

10 with the reported mode of solution. With this later computation, the Aiding Source 2 
may then apply the correction to the location reported by the SPS Receiver to obtain 
a more current location. This precise location then may be sent to the Call Taker 3 
in a Second Aiding Report 10. 

As an alternative embodiment of this aforementioned exchange, the SPS 

15 Receiver 1 reports the code phase differences to the Aiding Source 2. In this event, 
the Aiding Source 2 could compute a PV solution for the SPS Receiver 1 using a 
method like that of the SPS Receiver 1 . 

As described, the scheme requires the SPS Receiver 1 to compute 
pseudoranges without the benefit of having actual time synchronization data from the 

20 satellite signals. The use of this data is avoided because the code phases of the 
satellite signals are taken as ambiguous measurements of the differences in the 
times of transmissions of the satellite signals. This is accomplished by measuring 
the code phases at the same instant so that there is a common time of receipt. 
Moreover, the ambiguity resolution needed to convert code phase differences into 

25 pseudorange differences is achieved by utilizing the approximate location for the 

receiver obtained from the Aiding Source 2. When the approximate location is then 
combined with approximate time from a real time clock (accurate to say 1 minute) 
and current almanac data (less than about 2 months old in the case of GPS which 
provides range errors of less than about 30km), this permits approximate ranges of 

30 the satellites to be determined to better than half of a code epoch multiplied by the 
speed of light which is the ambiguity interval for the code phase differences. By 
combining the set of approximate ranges and the set of code phases, unambiguous 
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and precise pseudorange differences are derived without the need to synchronize 
the receiver to better than about one minute. 

Furthermore, consistent with the aiding data minimization objective, the SPS 
Receiver 1 is able to compute PV solutions without the need for being supplied with 
5 Doppler information. In traditional devices, such information is used to assist in 

restricting a carrier frequency search during signal acquisition or to predict changes 
in satellite range as an altemative to using the satellite ephemeris data. In one 
embodiment, the SPS Receiver 1 estimates Doppler information using stored current 
almanac data. 

10 Using this current almanac data, the SPS Receiver 1 estimates the Doppler 

frequencies of the satellites to an accuracy of better than about 250Hz for GPS given 
approximate locations of the satellites known to better than about one hundred 
kilometers. This is sufficient accuracy to achieve a rapid acquisition provided that 
the frequency offset of the reference oscillator of the SPS Receiver 1 is known to 

15 within a few Hz. To achieve this latter requirement, the reference frequency offset is 
estimated each time a PV solution is computed and thus it can be tracked. 
Moreover, to the extent that the reference frequency varies with temperature as well 
as aging and to the extent that a large change in temperature between PV solutions 
will lead to a degradation in acquisition performance, the SPS Receiver 1 , as 

20 described below, utilizes a method that copes with the change and ensures graceful 
degradation. 

To ensure the presence of current almanac data without the use of an Aiding 
Source 1 . the SPS Receiver 1 must be activated often enough and for long enough 
in the presence of strong signals to keep the data current. In the case of the GPS 

25 system, the almanac data must be less than 2 months old to remain current. To 

meet this goal under GPS. the SPS Receiver 1 must gather approximately 27 sets of 
orbital coefficients over a period of around 2 months. On average it would take 
around 20s plus signal acquisition time to gather one set and it would take around 60 
such gatherings to acquire all of the sets. Hence, if the SPS Receiver 1 was 

30 activated in the presence of strong signals approximately once per day on average 
for about 30s then the stored almanac would remain current. 
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Calculation of Almanac Aiding from Table of Orbit Coefficients 

It Is common practice to supply aiding information to a GPS receiver in the 
form of a full Ephemeris and its time of issue. One possibility is to have the receiver 
acquire its own Almanac data but Augmentation coefficients are supplied as aiding 
information. These coefficients permit the location and velocity computed from the 
Almanac data to be corrected to similar accuracy as provided by the Ephemeris 
coefficients. One disadvantage of this approach is that the receiver has to indicate 
the Issue Of Almanac to the aiding source so that the correct aiding data can be 
supplied. 

An alternative to this approach is to hardcode orbital coefficients for a number 
of the most common satellite orbits into a lookup table in the GPS receiver. The 
aiding source then supplies aiding to the receiver for a given number of visible 
satellites in the following form: 

Satellite ID 

Reference time of week (a required parameter of any orbit model) 
Look up number to a table of hard-coded orbital coefficients in the receiver 
Corrections to the orbital coefficients and/or additional orbital coefTicients 
and/or corrections to the output after using the orbital coefficients. 

Because of the compact nature of this data compared to full Ephemeris data 
or full Almanac data, it is feasible to structure the aiding in the form of broadcasts 
transmitted at regular intervals. It would be unnecessary for the receivers to request 
aiding since they would not need to supply the Issue Of Almanac and this is another 
reason why broadcast coefficients would be feasible. One advantage of using 
broadcasts is that this would significantly reduce system data traffic in the direction 
toward the network base center, call center or wireless web site. Another advantage 
of the scheme is that the receivers do not need to acquire Almanac data which 
means that they need not be operated in the open at regular intervals to do so. 

Calibration of a Local Oscillator 

In many Assisted GPS applications rapid acquisition is a key requirement 
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Reference oscillator drift is a major contributor to acquisition time as it causes the 
"Doppler" frequency search to be increased to allow for reference oscillator drift. By 
utilizing the precise signal framing of a digital communications link, it is possible to 
calibrate a local oscillator. This can be performed as follows: 

ko = Frequency of Local Oscillator in Hz to be corrected. 

fsF = Frequency of Signal Framing of digital communication link in Hz. 

ToIloc = lolerance of Local Oscillator before Calibration in PPM. 

TolgF = Tolerance of Signal Framing of digital communication link in PPM. 

Offset^o = Offset (in Hz) required to calibrate the Local Oscillator to tolerance 

ToIloc 

P = Period of time (in seconds) over which Local Oscillator clock cycles 

must be counted in order to calculate OffsetLo to Tol^cx:- 
NsF = Number of Signal Frames over which Local Oscillator clock cycles 
must be counted in order to calculate Offset^© to ToIloc- 
Nlo - Number of Local Oscillator clock cycles counted during P 
(or NsF Signal Frames). 

n = 1/n is the fraction of a Local Oscillator cycle able to be measured 
at the signal framing transitions. — 



The scheme involves counting local oscillator cycles and fractions thereof 
over a period precisely determined by a number of signal framing intervals. The 
measurement resolution relative to the local oscillator period is a 1/n fraction of the 
local oscillator period and is easily determined with the presence of any hardware 
clock whose frequency is approximately n times the local oscillator frequency f^o . 
As an example, the local oscillator may be fLo= 10 MHz and so the presence of a 50 
MHz clock with an accuracy of 1/2n Gust 10% in this case) or better will permit 
measurement of the local oscillator phase to the nearest 1/n of a cycle at the signal 
framing transition. The measurement is repeated at a later signal framing transition 
and the total number of local oscillator cycles and n^ of a cycle is determined. 
Assuming ToIsf « ToIloc . t^^©" 
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P = 10^ /(n X fLoxTolLoc) 
The minimum number of signal framing intervals to be used is given by 

5 Nsp = Roundup (fsp x P) 

where the "Roundup" function provides the next higher whole integer when (fgp x P) 
Is a non integer so that the correction meets or exceeds the accuracy required. 

10 The actual measurement period is I fsF 

where N^o cycles (accurate to 1/n of a cycle) are counted. 

The more accurate period for a single local oscillator cycle is Ngp/ (fgp x Nlo) 
The more accurate frequency is thus (fsp x Nlo) / Nsp 

15 

The frequency offset of the oscillator is then given by 

Offset^o = ^Lo - (fsF X Nlo) / Ngp 

20 This calibration offset can be used as a correction by the GPS receiver 

firmware when performing acquisition searches or it can be used to correct the 
oscillator frequency so as to minimize the offset. 

Weak Signal Acquisition/Tracking In Presence of Jamming Si gnals 

25 One of the problems associated with the use of weak BPS signals is that any 

SPS system has a limited dynamic range. In the case of the C/A Code signal of the 
GPS system, for example, any signal that is weaker by more than about 20dB than 
another signal that is also present may be jammed by the stronger signal. There are 
2 main effects of this jamming. First, while attempting to acquire the weaker signal 

30 using a sufficiently low threshold, the receiver's search sequence will be interrupted 
by frequent false alarms because the cross-correlations between the strong signal 
and the code generated in the receiver will be often be above the threshold. 
Second, although the receiver may be able to acquire and track the weaker signal, it 
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will be susceptible to large measurement errors caused by cross-correlation side 
lobes from the stronger signal dragging the tracking algorithm away from the true 
correlation main lobe. 

It is highly desirable to avoid the first of these problems and it Is essential to 
avoid the fatter problem as it can result in gross positioning errors of up to several 
kilometers. Thus, the SPS Receiver 1 should attempt to acquire strong signals 
before attempting to acquire weaker signals. In this regard, FIG, 2 outlines an 
example procedure for a GPS receiver that ensures that any strong signals are 
acquired first using a high threshold. A threshold for acquisition of any remaining 
needed signals will then be set, in the case of GPS, 20dB below the strongest signal 
acquired. While the FIG. 2 contains information pertinent to GPS system, its general 
application would be equally applicable to any other SPS system. 

Referring to FIG. 2, the device starts with the first request/response exchange 
in attempt to acquire the SPS Receiver's 1 approximate location from the Aiding 
Source 2. If the exchange is successful and the approximate location aiding data is 
received, the device sets its initial search parameters in step 1 1 with the goal of 
acquiring the strong signals. The parameters for the search are initially set to a high 
threshold and with a short integration period adequate for acquisition at the selected 
threshold. There is no prior knowledge of the code phases of the satellites and 
therefore the search is unrestricted. The reference frequency offset is assumed to 
be a prior value that was measured when the receiver was previously active (i.e. 
"old"). In subsequent step 12. with these parameters, the SPS Receiver 1 uses a 
multi-channel device to perform a concurrent search for the strong signals of ail 
visible satellites. 

In the GPS example, conducting such an unrestricted code search with 1 
channel per satellite at 1 chip per integration period using 2 arms per channel with 
0.5 chip spacing would take 1023/4 = 256 integration periods or approximately 1 
second. In most cases this would be adequate to acquire any strong signals 
present. 

However, In some Instances, the assumed frequency offset of the receiver's 
reference oscillator may have changed by a larger amount than can be 
accommodated for by the sampling rate of a correlator's output samples. Thus, if no 
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satellite signals are acquired during an attempt using the assumed reference 
frequency offset, then the offset would be adjusted in step 60 and the search 
continued without lowering the threshold. In this way, the offset would be changed 
systematically to effect a search over the possible frequency range. The systematic 
5 search would terminate on the acquisition of one or more satellites at the highest 
possible threshold. 

If no satellite signals were acquired over the possible frequency range then, in 
step 13, a lower threshold would be set (e.g. 6dB lower) and a longer integration 
period would be used (e.g. 4 times as long). The reference frequency offset would 

10 be set back to the previously assumed value and the frequency search would be 
restarted using these parameter values. 

In the original search to locate at least one satellite, if any strong signals are 
acquired and additional signals are still needed, then the measured carrier to noise 
ratio of the strongest signal acquired would be used to determine both the integration 

15 period to be used for a subsequent search in step 14 and the acquisition threshold to 
be applied during the search. For example, if a signal of greater than SOdBHz was 
acquired then the integration period for the subsequent search need not be any more 
than 32ms. This is because it is possible to detect signals of 30dBHz or more with 
an integration period of 32ms and the threshold would have to be set at 30dBHz or 

20 higher to avoid the dynamic range problems previously described. 

The reference frequency offset would also be estimated using the 
approximate location of the SPS Receiver 1 together with the almanac data and the 
measured carrier frequency of the acquired signal. This should obviate the need for 
further frequency searches. 

25 Since at least one satellite has been acquired, a restricted search regime can 

be conducted for the remaining satellites using one channel per satellite since the 
approximate code phase differences between all of the remaining visible satellite 
signals and the first one can be estimated. If the accuracy of the approximate 
position estimate is assumed to be ±10km then the code phase differences can be 

30 estimated to within ±25, approximately. In the GPS case this would permit all of the 
remaining satellites to be acquired within 50 * 0.128s or 6.4s. However, the search 
could be terminated once sufficient satellites had been acquired to permit the 
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location to be determined. 

In the GPS example, assuming that the second search did not fall, the 
maximum time taken to acquire enough satellites would thus be 1s for the first 
search for strong signals plus 4s for the second search for the first satellite plus 6.4s 
5 for the subsequent search for the remaining satellites. This adds up to 1 1 -4s. 

However, typically, acquisition would take less than this (e.g. 1s + 4s + 3.2s or 8.2s). 

Cancellation of Cross Correlations 

A preferred version of the basic idea of canceling strong signals deals with the 

10 problem that the signals are represented with very low precision at the input to a 
correlator and, hence, any scaling of the signal can only be extremely crudely 
performed. A preferred variation which allows the scaling to be performed at a point 
where the signal is represented by a higher number of bit samples, e.g. 10 bit 
samples and scaling is more feasible is the following: 

1 5 The concept involves hardware provision in the correlator for cancellation. 

Fig. 7 is a schematic diagram showing how a canceler may be integrated into a 
correlator Each channel or certain channels may be reconfigured to provide a 
canceling service for others. Alternatively, these canceling elements can be 
dedicated for the purpose. For each finger of each channel for which canceling is to 

20 be performed, a canceling finger is required. 

The canceling finger takes one input from the code generator input to the 
prompt finger of a channel tracking a jamming signal. It takes the second input from 
the code generator input of the finger for which cancellation is to be performed. By 
correlating these two signals together, the canceling finger generates an estimate of 

25 the cross-correlation between the two C/A-Codes involved. 

This is then scaled in proportion to the magnitude of the jamming signal which 
is obtained from the correlation achieved In the channel tracking the jammer The 
result is an estimate of the value of the cross-correlation experienced in the channel 
for which cancellation is required. Cancellation is then achieved by simply 

30 subtracting the estimated cross-correlation value from the correlation measured. 

This embodiment has the advantage that the scaling and cancellation are 
performed at a point where the signal is represented by approximately 10 bits. For 
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this to be so the bandwidth of the signal must have been reduced by a factor of 
around 2^ or 256 times through a process of decimation (integration and down- 
sampling). It is feasible if all channels have been downconverted identically and 
therefore the downconversion and code mixing steps can be as shown in Fig. 7. The 
last channel specific downconversion step is small but large enough to account for 
the Doppler variation of plus or minus 4KHz. 

Simultaneous Code Phase Determination of Multiple Satellite Signals 

To achieve the aiding scheme as described, an SPS Receiver 1 of the 
present invention requires the ability to determine the code phases for multiple 
satellite signals at the same instant. In this regard, FIG. 3 depicts one such SPS 
Receiver 1 . Generally, the SPS Receiver 1 can be broken down into roughly three 
parts. The device has a front-end circuit 17, three or more correlators 18 and a 
microprocessor 20 with memory. The following describes the functions of each. 

Generally, the front-end circuit 17 serves as the initial signal processor as 
follows. The front-end circuit 17 amplifies, filters, down converts and digitizes the 
signal from an antenna so that it is suitable for processing in a digital correlator 18 
and such that the signal to noise and signal to interference ratios are minimized 
subject to economic and practical realization requirements. The front-end output 19 
of the front-end circuit could be a complex signal centered at tens of KHz (in the 
case of GPS) or a real signal centered at around 1.3MHz or higher. The sampling 
rate would typically be several MHz and the digitization would be at least 2 bits per 
sample. In the preferred embodiment an AGC circuit keeps the level of the digitized 
signal constant. Since the true signals are spread over 2MHz in the case of GPS 
and are weak signals in any case, this signal is dominated by noise and the AGC 
maintains a constant noise level at the output of the front-end. 

Hardware conrelators 18 each representing a processing channel for a 
particular satellite signal are used separately to further process the front end output 
19 under the control of the microprocessor 20. Within each correlator 18 a further 
down conversion 21 (quadrature in this case) to neariy DC is performed based on 
the estimated Doppler offset of a particular satellite signal and the estimated offset of 
the crystal oscillator reference frequency driving the correlator. Then, the resulting 
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complex down converted signal 22 is mixed with (i.e. multiplied by) a real binary 
pseudorandom code signal 23 chosen to match that of a particular satellite signal 
and generated by a code generator 24. The code generator 24, controlled by the 
microprocessor 20, generates the pseudorandom code signal 23 at a selected rate 
5 set to match the estimated signal Doppler offset given the estimated crystal oscillator 
offset. 

The code generator 24 also generates a late pseudorandom code signal 25 
that is the same as pseudorandom code signal 23 but at a fixed lag with respect to 
the former. This late pseudorandom code signal 25 is also mixed with the down 

10 converted signal 22. The resulting mixed signals 26 are then separately processed 
by decimators 28. Decimators 28 low pass anti-alias filter and down sample the 
mixed signals 26 to a reduced sampling rate. In the case of GPS, the reduced 
sampling rate is approximately IKHz. This sampling rate may be derived from the 
local code rate such that a single sample will be obtained for each code epoch. 

1 5 However, this is not essential. 

When searching in code, the processor 20 either causes the code generator 
24 to step instantaneously by the required amount at the start of each integration 
period or changes the code frequency by a known amount for a precise period of 
time so as to effect a rapid step in the code lag. This is the preferred embodiment 

20 although in an altemative scheme the code frequency may be deliberately offset 
while searching so that the code slews continuously relative to that of the incoming 
signal. 

When tracking a satellite signal, in this embodiment, the processor 20 
constantly adjusts the code lag as just described so as to keep the pseudorandom 

25 code signal 23 and late pseudorandom code signal 25 from the code generator 24 
running one ahead (eariy) and one behind (late) the code of the incoming signal. In 
other embodiments, the code generator 24 may generate a third signal (prompt) that 
is kept running synchronously with the code of the incoming signal or, indeed, there 
may be several more signals spanning a lag interval of up to 1 chip (smallest code 

30 elements) early and late of the incoming code. 

The correlator output samples 29 are read into the processor 20 where they 
are further processed by a signal processing algorithm described later In this 
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specification to estimate the amplitude, frequency and phase of the carrier signal. 
Then. If data is to be extracted because the signal is strong enough for that, the 
phase and frequency are utilized by a separate algorithm that operates on the raw 
samples to extract the data. Methods for the extraction of this data will be obvious to 
one skilled in the art. 

The frequency of the tracked carrier signals are then used to estimate the 
Doppler offset of the carriers and the crystal oscillator offset. The former Doppler 
offset values are subsequently used to estimate the velocity of the receiver (and the 
vehicle In which it may be traveling). 

The amplitude of the eariy and late correlator output samples 29 represent 
estimates of the carrier to noise ratio of the satellite signal since the noise level is 
maintained constant by the AGC of the front-end. When performing a satellite signal 
search, the amplitudes are compared to a threshold to determine If the signal has 
been detected. If It has, then an acquisition procedure is commenced. The steps of 
an appropriate acquisition procedure will be obvious to one skilled in the art. 

When tracking, the code phase is adjusted as discussed eariier in order to 
keep the average amplitudes of the correlator output samples 29 equal to each 
other. A similar, but more complex algorithm can be applied If 3 or more correlator 
output samples 29 are present. The nature of these control algorithms will be 
obvious to one skilled in the art. 

At the end of each integration period the code phase 30 for each correlator 18 
are simultaneously latched by a latch element 31 within the hardware correlator. 
The resulting signal represents the code phase measurement 32, These code 
phase measurements 32 are then made available to the processor 20. The . 
processor 20 then applies a smoothing algorithm to the code phase measurements 
32 together with the carrier frequency estimates. This algorithm Is used to reduce 
the random error in the code phase measurements 32 over time by making use of 
the precision in the carrier frequency estimates to predict the changes in code phase 
from one integration period to the next. The algorithm also filters the carrier 
frequency estimates to reduce their random errors. The carrier smoothing algorithm 
is described later in this document. 

After several seconds of smoothing, (five seconds in the preferred 
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embodiment), the carrier smoothed code phase measurements and the filtered 
carrier frequency estimates are passed to the location solver which estimates the 
SPS Receiver 's 1 location and velocity. The algorithm makes use of precise 
satellite position data and the approximate location received from the aiding source. 
5 This calculation is described in more detail later In this document The signal 

processing, carrier smoothing and location solving algorithms are all executed by the 
processor 20. 

Modified Search Engine 

10 Search engines vary from chip to chip. One well known type (Solution 1 ) 

involves 12 channels with 20 fingers per channel and these channels can be ganged 
to provide 240 fingers. In this ganged mode each finger has 8 DFT bins thereby 
providing the equivalent of 1920 fingers all searching for the one satellite. Another 
well known type (Solution 2) Involves a search engine with approximately 16,000 

15 fingers which can be used to search for up to 8 satellites simultaneously across the 
full code space but with only one frequency bin. Other types Incorporate limited 
search engine capacity together with correlators. 

in ail cases the search engine hardware supports coherent Integration, usually with 
variable. Integration period, followed by non-coherent Integration over a variable 

20 period. Fig. 8 illustrates a preferred 2 bin search engine architecture. The 
coherent interval determines the binwidth but It also determines the coherent 
processing gain. The squaring loss incurred at the end of the non-coherent 
integration step depends on the signal-to-noise ratio prior to squaring and therefore it 
is desirable to maximize the coherent interval so as to minimize the squaring loss. 

25 However, it is necessary to obtain sufficient bandwidth to cope with residual TCXO 
error and Doppler error due to Initial positioning error. With longer coherent 
integration periods this may require frequency searching unless multiple frequency 
bins are employed. 

In a wireless location application a minimum of around 200Hz bandwidth Is 

30 required to allow for up to 20km cell radius and .1PPM residual TCXO error after 
network compensation (if employed). 

Since the Solution 2 hardware does not incorporate a correlator it is 
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necessary for all of the satellite search engines to commence simultaneously at the 
same codephase so that the codephase measurements are meaningful in a relative 
sense. This rules out sequential frequency searches and thus limits the coherent 
interval that may be employed to 5 ms so as to provide 200Hz bandwidth. If only 4 
5 satellites are acquired simultaneously then 2 bins can be produced for each and the 
coherent inten^al can be increased to 10s. However, with the fading typical of weak 
signal environments, this would not produce good performance. 

Where correlators are employed, tracking firmware may be used to track 
satellite signals following acquisition using the search engine. This provides more 

10 freedom in the optimization of the search algorithm. Longer coherent integration 

periods may be employed because the search may be undertaken sequentially and 
satellites may be acquired sequentially. Once the first satellite has been acquired 
the search range may be restricted in code space to reduce acquisition time or the 
required search engine capacity. 

15 For an architecture like Solution 1 the Search Engine mode is employed for 

the first satellite after which it is far quicker to use the correlators with acquisition 
firmware based on the algorithms disclosed in this document to perform in parallel a 
restricted sequential code search for all of the remaining satellites. The search may 
be repetitive to allow for fading. Several of the 20 fingers may be used in parallel to 

20 speed the search process although the processing load may be prohibitive if more 
than 6 fingers are used. (The fingers are used in pairs to minimize processing load 
and to optimize detection performance.) 

Searching using conventional correlators with the acquisition algorithms of the 
invention is a sequential process. However, because the integration periods 

25 required for given sensitivity are far less using the processing scheme of the 

invention than using a search engine, the choice of correlators or Search Engine for 
acquisition depends on the available search engine capacity. For the first satellite 
the search engine is invariably quicker. For subsequent satellites it may be quicker 
to use correlators. 

30 A special search engine configuration may be employed to integrate with the 

processing algorithms of the invention. This allows the search engine to take 
advantage of the greatly reduced integration periods of the invention at the cost of 
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more frequency bins. The fingers are used in groups of 30 to provide consecutive 
DFT bins to allow for the data bandwidth and frequency uncertainty with the 
narrower binwidth corresponding to a coherent integration period of 128ms. Prior to 
magnitude squaring, an extra processing block is required but the extra processing 
5 capacity required represents an additional overhead of only 10% of the coherent 
processing capacity requirement. After non-coherent integration the peak bin and 
adjacent bins are passed to the firmware for additional processing to estimate carrier 
magnitude and frequency for detection and handover to the correlator. Fig. 9 
illustrates the modified search engine architecture. 
10 Using 1ms minimum and 19ms maximum coherent integration periods, the 

non-coherent integration time required for detection of a signal using a conventional 
search engine is given as in Table 1. (A choice of 19ms provides a long coherent 
integration period below 1 bit and ensures that the bit periods slide continuously 
through the coherent integration periods. This results in a loss of 1.2dB but avoids 



15 any requirement for data synchronization.) For 5ms coherent Integration the non- 
coherent Integration periods required would be much longer. 
TABLE 1 - SEARCH ENGINE INTEGRATION TIME VS CONDITIONS 





CONDJTION 


RF SIGNAL 


IFC/No 


INTEGRATION 






LEVEL (Dbw) 


(dBHz) 


PERIOD (ms) 




Open skies, directly above 


-150 


54 


1 


20 


Open skies, high elevation 


-155 


49 


1 




Open skies, mid elevation 


-160 


44 


2 




Open skies, low elevation / 


-165 


39 


6 




indoors, very strong signal 










Open skies very low elevation 


-170 


34 


19 


25 


/ indoors, strong signal 










Indoors, moderately strong 


-175 


29 


88 




signal 










Indoors, weak signal 


-180 


24 


599 




Indoors, very weak signal 


-185 


19 


5270 


30 


Indoors, extremely weak 
signal 


-190 


14 


51500 



The lower the signal to noise ratio prior to squaring the higher the squaring 
loss and this is the reason that the required integration period rises exponentially. 
This, in turn, is the reason why it is essential to acquire satellites in parallel and 
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hence why multiple search engines are needed. 

However, the integration period required for the processing scheme of the 
invention is less because much longer coherent integration periods are employed 
thereby reducing the squaring losses. 
5 Table 2 compares the integration periods and acquisition times of 3 schemes 

for very weak to extremely weak signals. The hardware employed in the 3 schemes 
is a search engine of around 16,000 fingers plus, say. 12 correlators. This is just an 
example. Similar analyses can be performed for other search engine configurations 
and the results in all cases favor the use of the algorithms of the invention either 
10 running in firmware using the correlators or in the form of a hardware mechanization 
as in Figure 9. 

The first scheme uses the Search Engine only with 5ms coherent integration 
periods. The second scheme uses the search engine sequentially with 19ms 
integration periods for acquisition and the correlators for tracking. Only 2 steps are 

15 required as all remaining satellites can be acquired in parallel after the first satellite is 
acquired. The full search engine is used to search for 2 satellites in parallel in the 
first step. The third scheme uses a DFT search engine and a hardware 
mechanization of the algorithms of the invention for sequential acquisition along with 
correlators and firmware based on the algorithms of the invention for tracking. Three 

20 search steps are required, 2 for the first satellite and one for the remainder. 

TABLE 2 - INTEGRATION PERIOD & ACQUISITION TIME FOR SEARCH ENGINE 
Vs THE INVENTION 



Sensitivity 


SCHEME 1 


Search 


SCHEME 2 


Worst Case 


SCHEME 3 


(dBW) 


Search 


Engine 


Acquisition 


Integration 


Modified 




Engine 


integration 


Time With 


Period (s) 


Search 




Acquisition 


Period With 


19ms 


For The 


Engine 




Time With 


19ms 


Coherent IP 


Invention 


Acquisition 




5ms 


Coherent IP 


(s) and 




Time (s) 




Coherent IP 


(s) 


Tracking 








(s) 










-185 to - 


38.2 


15.2 


30.4 


4.2 


12.6 


187 












-186 to - 


60.3 


24.0 


48 


5.4 


16.2 
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188 












-187 to - 
189 


95.1 


37.8 


75.6 


8.3 


24.9 


-188 to - 
190 


154.7 


60.2 


120.4 


13.1 


39.3 



To allow for 20km of initial position uncertainty we require a search range of at 
least 60 chips after acquisition of the first satellite. This deternnines the number of 
search engine fingers required after acquiring the first satellite. This is 480 per 

10 satellite for scheme 2 and 3,600 for scheme 3 which means that, in both cases, four 
or more of the remaining satellites can be acquired in parallel. 

As we have seen, Scheme 1 requires the longest acquisition time and hence 
the highest energy requirement for the first fix. Moreover, subsequent fixes require 
equally as much energy as the satellites must be acquired from scratch again. Using 

15 the invention for tracking, however, the energy required for subsequent fixes can be 
greatly reduced in 2 ways. Firstly, the power hungry search engine can be turned 
off. Secondly, the entire receiver can be duty cycled over quite long periods to 
achieve considierable savings. A Timer must be used to turn the receiver back on 
after a period precisely known to within a few microseconds. Then the code 

20 generators can be shifted close to the tracking codephase prior to a brief 

reacquisition search using the invention implemented in firmware. Since only the 
correlators are used the energy requirement is substantially less than would be the 
case for the first fix. 

25 Satellite Signal Processinq Aiaorithms 

As previously noted, special processing of the weak satellite signal is 
required in order to determine its code phase differences for purposes of calculating 
a PV solution from them. To this end, the invention utilizes an algorithm to measure 
amplitude, frequency and phase of the satellite signal. The invention also applies a 
30 smoothing procedure to the code phase measurements and the carrier frequency 

estimates to reduce random errors occurring over time. Finally, a formula is applied 
to convert the code phase differences into a precise position and velocity solution. 
Each algorithm is addressed in turn. 
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(1 ) Estimating Amplitude. Frequency and Phase of Weak Satellite Signal 

Any SPS Receiver 1 employing a hardware correlator 18 with several signals 
in each channel is required to estimate the amplitude of the carrier in each of those 
several signals in the presence of the data. In a conventional receiver a phase 
5 locked loop or a frequency locked loop and a delay locked loop control the frequency 
of the final down converter 21 and the code generator 24 respectively. In the weak 
signal case however, the signals are too weak to permit lock-in without the use of 
some sort of aiding from an auxiliary algorithm. The signal processing algorithms of 
the present invention could be used as auxiliary algorithms for acquisition or they 

10 may be used independently of any phase-locked or frequency-locked loop for both 
acquisition and tracking. 

With strong signals it is possible to detect a signal from individual correlator 
output samples 29 with sufficient reliability to permit lock-in. In the weak signal case, 
the correlator output samples 29 are so noisy that the signal is indistinguishable from 

1 5 the noise unless the correlator output sampling rate is extremely low (e.g. 8 Hz for 
GPS). However, to avoid the Decimators 28 filtering out the residual carrier at such 
a low sampling rate, the Doppler frequency and the crystal oscillator offset would 
have to be known to an impractical high precision. Accordingly, a higher correlator 
output sampling rate needs to be retained and a suitable algorithm is required to 

20 estimate the amplitude of the residual carrier signals in the code phase 
measurements 32. 

FIG. 4 depicts a flow chart for an algorithm used to measure amplitude, 
frequency and phase of the signal from each of the eariy and late mixed signals 29 
of the correlator 18 according to one embodiment of this invention. Referring to FIG. 

25 4. the procedure involves the use of a Fast Fourier Transform 33 applied to a block 
of samples from the code phase measurements 32 of the correlator 18. The effect 
of this algorithm is to compress the residual carrier signal into a few bins of the FFT 
output 34. Thus, whereas the satellite signal is undetectable when its energy is 
spread across all of the time domain samples contained in the code phase 

30 measurements 32, it is detectable in the bins of the FFT output 34. However, its 
amplitude is not readily estimated because the data modulation splits the satellite 
signal between several adjacent bins in a semi-random manner depending on where 

30 
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the transition falls in relation to both the integration period and the phase of the 
residual carrier at that Instant. 

Nevertheless, having detected the probable presence of a satellite signal it is 
possible to reject much of the noise by applying a window operation 35 on the bins of 
5 the FFT output 34 centered on the peak value. More distant bins may be discarded 
completely. This is simply a filtering operation and significantly improves the signal- 
to-noise ratio prior to non-linear processing 36 of the window*fiitered signal 37 to 
eliminate the data transitions. The remaining windowed bins in the window-filtered 
signal 37 may be processed in one of several ways in order to estimate the 

10 amplitude in the presence of the noise and data as follows: 

1 . The discarded bins may be zero filled and the complete set of bins can then 
be inverse transformed back to the time domain. This would result in a set of time 
domain samples with significantly improved signal to noise ratio which may then be 
processed such as to effect iock-in of a phase locked loop or delay locked loop. 

15 2. An inverse transform may be effected as described above and the time 

domain samples may be squared to remove the data. The amplitude of the resulting 
signal could then be estimated by transforming back into the frequency domain using 
a DFT to obtain a few bins centered on the previous FFT peak and applying any 
estimation algorithm suitable for estimating the amplitude and frequency of a cissoid 

20 embedded in noise based on the values of several FFT bins. It is important to note 
that one effect of squaring the signal will be to double the frequency of the residual 
carrier signal such that it may be aliased by the samples. This ambiguity will need to 
be resolved when determining the frequency. 

3. The vector of remaining windowed bins may be autoconvolved to remove 

25 the data in an autoconvolution process 36. This is equivalent to squaring in the time 
domain but, if the number of bins in window filtered signal 37 is small compared to 
the size of the FFT 33, then the autoconvolution process 36 may be less 
computationally costly than the process described in option 2 above. The 
autoconvolved samples 38 can then be processed by any estimation algorithm 39 

30 suitable for estimating the amplitude and frequency of a cissoid embedded in noise 
based on the values of several FFT bins. Again it is important to realize that the 
frequency corresponding to each bin has been doubled by the autoconvolution 
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process 36 and the bin width has thus been effectively halved for the purposes of 
frequency estimation. 

An embodiment of the invention employs option 3 for estimation of the 
amplitude and residual carrier frequency. The RF carrier frequency is estimated as 
follows: 

Fc = Fd1 + Fd2 + (Np1 + (Np2 - Nnom + )/2)*Fbin - 4^ Fxo*Fnom/Fxo 
Where: 

Fc is the RF carrier frequency; 

Fd1 is the total frequency shift due to the front-end down conversion; 

Fd2 is the frequency shift due to the down conversion in the correlator; 

Np1 is the bin number (between -N/2 and (N/2+1 )) for an N-point FFT) of the 
peak in the FFT (which is the center bin of those extracted for performing the 
autoconvolution); 

Np2 is the bin number of the peak bin in the autoconvolved bins computed 
(which is the center bin of those extracted for estimating the amplitude and 
frequency); 

Nnom is the bin number of the nominal peak bin in the autoconvolved bins 
computed (i.e. corresponding to zero lag); 

N is the frequency adjustment (in bins and fractions of bins relative to Np2) 
estimated by analysis of several adjacent bins of the autoconvolved bins computed; 

Fbin is the bin width of the original FFT; 

Fxo is the crystal oscillator offset from its nominal frequency; 

Fnom is the nominal RF carrier frequency of the signal; and 

Fxo is the nominal crystal oscillator frequency. 

While the signal is being tracked it is possible to reduce the computational 
load of the signal processing, if this is desirable, by utilizing the precise value of the 
residual carrier frequency as estimated in the previous integration period. For 
example, only those bins used in the previous integration period need be computed. 

Ensemble Averaging Scheme 

The coherent algorithm is limited in sensitivity owing to the fact that gains 
resulting from the reduction of binwidth are offset by squaring losses. A more refined 
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averaging scheme has been tested and seems preferable. 

The basic concept of this scheme is to extract the signal power spectrum from 
the noise by averaging. Because the phase information is lost in computing the 
power spectrum, the MacLeod estimation algorithm must be replaced by the less 
5 efficient (in an estimation performance sense) Jain algorithm for estimating 
amplitude and frequency for detection and carrier tracking. However, the 
consequent degradation in estimation performance is more than offset by the 
improved sensitivity possible. Nevertheless, when shorter integration periods are 
required, the MacLeod scheme should be retained. 
10 To avoid excessive memory usage for storing FFT arrays, the scheme must 

operate in the fashion of a moving average rather than by block averaging. This 
implies running a filter on the squared magnitude of the autoconvolution array for 
each channel. For smooth startup, such a filter should be of a simplified Kalman 
style. 

15 The first difficulty to be faced in implementing such a scheme is that the 

autoconvoluttons are normally performed on windowed FFT values and the window 
is centered on the peak value of the FFT during acquisition. This is not the case 
when tracking as the center of the window Is predetermined. Hence the scheme 
outlined above is very suitable for tracking but needs augmenting for acquisition as 

20 the peak picking of the individual FFTs will not be reliable at the low signal levels to 
be targeted. 

The solution proposed for this problem is to filter the squared magnitudes of 
the FFT bins in the same fashion as for the autocon volutions when acquiring. One 
way to facilitate this would be to dwell for several FFT periods at each code search 

25 step. An alternative approach would be to step very slowly at each FFT period, (ie 
for N-fold averaging apply 1/Nth of the normal search step size.) The 
autoconvolution and later processing could be performed at each step as normal 
with no ill effects. Since the steps will be very small, the signal spectra will be 
present for many steps, consistent with the averaging period, as required. There will 

30 be a lag between the signal appearing and actually being incorporated into the 
autoconvolution average because the FFT average may not produce the correct 
peaks initially. However, this should not pose a serious difficulty and can be taken 
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account of by adjusting the capture process as necessary. 

Finally, we need to consider the Impact of this new scheme on the tracking 
control algorithms. The approach proposed is to perform partial control adjustments 
consistent with the averaging period. For example, if the averaging period is N FFTs 
5 then the feedback at each should be applied with 1/Nth of the normal gain. 

Simulations of block averaging were conducted to confirm the increases in 
sensitivity expected. An FFT size of 128 was used with windows of 15 bins before 
and after autoconvolution. The Jain algorithm was used to estimate the amplitude. 
Statistics were accumulated over 1,000 runs. 

10 The mean, standard deviation and 99th percentile value of the amplitude 

estimates were computed with no signal present. Then the signal amplitude required 
to make the tenth percentile value equal to the No Signal 99*^ percentile was 
determined by iteration. The mean and standard deviation of the amplitude estimate 
were recorded for this actual signal amplitude. This whole process was repeated for 

1 5 many different averaging periods. 

Table 3 shows the results obtained form the simulations. The predicted 
sensitivity is consistently neariy 3dB better than that predicted using a detailed 
spreadsheet analysis. This Is because the simulation does not simulate the losses 
due to quantization and correlation offset which amount to 1 .8dB before squaring 

20 and more after squaring. After taking this into account the simulations confirm the 
analysis well. 

One problem that is evident is that the signal amplitude estimates drift further 
and further away from the true amplitudes. At the same time the noise amplitude 
standard deviation decreases more rapidly than the mean noise amplitude. This 

25 suggests that the amplitude estimates could all be adjusted by subtracting a factor 
that depends on the value of N, the averaging length. The table shows the effects of 
applying such an adjustment based on the empirical results. This will not affect 
delectability and has no other adverse characteristics. Therefore, it is proposed to 
use a look up table based on these results for the purpose of applying the 

30 adjustment (Note that the adjustment for N=1 should be 0.0.) 
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The idea is to filter the autoconvolution and the FFT values. This means that 
there has to be a stored array for each channel containing the filtered value. One 
may as well use the same structure as for AutoConv except that the .Bins are real 
rather than complex. Use a simplified Kalman style approach to starting the filter 
5 but stop the convergence when the required Gain factor ( <1 ) is reached. 

The required Gain factor would be derived from a configuration parameter. It 
corresponds to the required sensitivity improvement in terms of signal power. If it is 
set to 1 then no filtering is needed. 

Note that, when searching, the step size at each step should be a fraction 
10 (also 1/G) of the normal chip step rather than just one integration period. 

It is proposed to use the required sensitivity in dBW as the configuration parameter. 
We could have one for tracking and another for acquisition. The Table 4 shows how 
G (= 1/Gain) relates to the required sensitivity. 



1 5 TABLE 4 - FILTER GAIN VERSUS SENSITIVITY 

Required Sensitivity (dBW) Required Sensitivity (dBHz) G 

-182 22 1 

-185 19 5 

-186 18 7 

20 -187 17 14 

-188 16 19 



One wouldn't choose -188 for acquisition using a conventional correlator because of 
how slow the acquisition would be. However, one could choose it for tracking based 
. 25 on the theory that, if you need such sensitivity then you must be indoors and moving 
at walking pace. Alternatively, a modified search engine, as discussed elsewhere in 
this document, may be used for acquisition using this scheme. 

Finally, G should be limited by the 20dB threshold limitation also. This should 
override the configuration parameter. 
30 The following pseudocode shows how the filter would be applied for the 

autoconvolution. 

if (Tracking) 
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Greq = Gtracking 

else 

Greq = Gacq; 
jf (Greq < 1 ) 
{ 

if ( first time ) 
{ 

FiitAutoConv(Clian).Bins = AutoConv.Bins 
G = 1 

} 

else 

{ 

K = 1/(1+G) 

FiltAutoConv(Chan).Bins = FiitAutoConv(Chan).Bins ->- 
K*(Magsquared(AutoConv.Bins) - FiltAutoConv) 

if (Greq <G) 

G = (1-K)*G 

} 

FiltAutoConv(Chan).Peak = FindPeak(FiltAutoConv(Chan)); 

FiltAutoConv(Chari).FFTPeak = FFTPeakA.Peak; 
/* */ 

/* Perform Jain Signal Analysis. 

Jain(&FiltAutoConv(Chan),pEstSignal); 

} 

else 

{ 

The algorithmic enhancement to the signal processing scheme avoids the 
need for impractical FFT lengths without unnecessarily limiting the sensitivity. The 
technique involves averaging of the power spectrum at both the FFT and 
autoconvolution stages. To keep memory requirements to practical levels, the 
scheme can be implemented using a filtering approach. 
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A simulation study has confirmed the effectiveness of the scheme in terms of 
improving the sensitivity with performance close to that predicted by detailed 
spreadsheet calculations. The study also yielded data that has allowed an amplitude 
adjustment calculation to be devised involving the use of a lookup table. 
5 An implementation scheme has been proposed and implemented in prototype 

form with a significant improvement in sensitivity resulting. 
(2) Signal Smoothing 

While the signal is tracked the code phase and frequency estimates may be 
refined by a carrier smoothing process prior to being passed to the location solver. 

10 FIG. 5 depicts the flow of such a process. The algorithm illustrated is a block 
computation applied to the stored estimates from several integration periods 
(arbitrarily labeled from J to N) to compute a single set of measurements 
representing refined measurements for the final integration period. 

In the preferred embodiment, the algorithm processes differences between 

15 the estimates from one satellite and those from all the others rather than processing 
the absolute estimates for individual satellites. The reason for this is that in the 
preferred embodiment of the location solver, differences rather than absolute 
estimates are used because differences can be processed without the need to take 
account of the frequency offset of the reference oscillator. 

20 All of the frequency difference estimates 40 corresponding to a particular 

satellite signal are simply averaged in step 42. The average carrier frequency 
difference 43 is then used to predict forward in prediction step 44. This prediction 
step 44 is used for all but the latest code phase differences 41 for that satellite to the 
latest measurement instant before averaging them in step 45. The prediction step 

25 44 uses an estimated Doppler offset of the code 47 determined from step 46 which is 
based upon the difference between the Doppler offsets of the carrier frequency (Fcl- 
Fc2). With the estimated Doppler offset of the code 47. step 44 is used to predict 
forward the code phase differences 41 by N integration periods using the following 
formula: 

30 ^p1 - ^p2 = Frac( ml - m2 + N*Tip/(Te* Fnom /( Fc1-Fc2)) 

where: 

pi - 4^ p2 is the predicted code phase difference 48 expressed as a fraction 
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of a code epoch; 

the Frac function yields the fractional component of the real argument; 
4^ m1- ^^m2 is the measured code phase (41) to be predicted fonvard; 
Tip IS the nominal integration period; 
5 Te is the code epoch period as determined from the tracking algorithm. 

All other quantities are as previously defined. 

The predicted code phase differences are then averaged in step 45 to 
produce an average code phase difference 49. 

(3) Position and Velocity Solution Computation using Code Phase Differences 
10 As previously discussed, a PV solution is computed by the SPS Receiver 1 

after code phase differences have been processed by the aforementioned smoothing 
algorithm. FIG. 6 illustrates the PV solution process used by a preferred location 
solver. The location solver, which computes the location and velocity using the 
refined average carrier frequency differences 43, the average code phase 
15 differences 49, the precise satellite positions 52 and the approximate receiver 
location 54 as follows. 

In step 50, approximate ranges 51 to all of the satellites for which there are 
measurements are computed. Since the satellite positions 52 are supplied by the 
Aiding Source 2, this step simply Involves the vector difference of the Cartesian 
20 coordinates of the satellite positions and the approximate location 54, which was 
also supplied by the Aiding Source 2. The vector magnitudes of these vector 
differences are approximate ranges 51 . 

In step 55, the epoch ambiguities are resolved and the average code phase 
differences 49 are converted into pseudorange differences 56. All of the 
25 approximate ranges 51 are subtracted from the approximate range of a reference 
satellite (as selected for computing the code phase differences) to give approximate 
range differences. These are saved for later use as well as being used in resolving 
the epoch ambiguities according to the following formula: 

P1-P2 = int[(R1-R2)/c*Te-( Q 1- Q 2)+0.51 + ( Q 1- Q 2) 
30 Where: 

C IS the speed of light; 

Te is the nominal epoch period; 
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P1 and P2 are the pseudoranges 56 of satellites 1 and 2; 
R1 and R2 are the approximate range estimates for the same two satellites; 
Q 1- Q 2 is the code phase difference between the same two satellites. 
In step 58, the range rate differences 57 are calculated from the Doppler 
5 affected average carrier frequency differences 43. The range rate differences 57 are 
computed according to the formula: 

d(R1-R2)/dt = -c*(Fc1-Fc2)/Fnom 
where: 

c is the speed of light and 
10 d(R1-R2)/dt is range rate difference 57 between satellite 1 and satellite 2. 

In one aitemative embodiment of the invention, the current location estimate 
60 and current velocity estimate 61 are computed in step 59 by a method similar to 
that for a single update of a Kalman navigation filter in a traditional SPS receiver. In 
fact, If time allows, a true navigation filter can be run for several updates to further 
1 5 refine the estimated PV solution. (However, in this case the satellite positions will 
need to be updated from the Aiding Source 2.) The Kalman gain matrix, K, is 
given by the well-known equation: 

K = PM^(MPM^+R)'' 

and 

20 X = X,^,j + K(Y - Yppeo) 

where Y is the measurement vector and X Is the solution state vector 
containing the location and velocity estimates. 

Using this equation, the initial state vector, X|n,t, is set to the approximate 
location from the aiding source with zero velocity. The first prediction vector, Yp^Eo. 
25 is set to the approximate range differences derived from the approximate ranges 51 
and zero for the range rate differences 57. 

The state covariance matrix, P, is initialized to a diagonal matrix with the 
entries representing the estimated variances of the approximate location and velocity 
estimates. The location variance estimates may be obtained from the Aiding Source 
30 2 or a fixed value could be used. The initial velocity estimate is zero and its variance 
depends on the application. The measurement variances can be estimated as In a 
conventional receiver except that the measurement variance matrix, R, is no longer 
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diagonal. The fact that differences with a reference satellite were computed means 
that the covariance terms t)etween the satellites are half of the value of the variance 
estimates rather than zero as in a conventional receiver. The approximate location 
(54) and the satellite positions (52) can be used to determine the direction cosines 
5 and the direction cosine differences form the rows of the measurement matrix, M. 
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